Magnetic Frustration on a Kagome Lattice in 
RsGasSiOM Langasites with R = Nd, Pr 



P Bordet^ I Gelard^ K Marty \ A IbanezS J Robert^, V 
Simonet^, B Canals^, R Ballou^, P Lejay^ 

1 Laboratoire de Cristallographie, CNRS, B.P. 166, 38042 Grenoble Cedex 9, France 

2 Laboratoire Louis Neel, CNRS, B.P. 166, 38042 Grenoble Cedex 9, France 

^ Centre de Recherches des Tres Basses Temperatures, CNRS, B.P. 166, 38 042 
Grenoble Cedex 9, France 

E-mail: pierre . bordetOgrenoble . cnrs . f r , ballouSgrenoble . cnrs . f r 

Abstract. In the R3Ga5SiOi4 compounds, the network R of rare earth cations form 
well separated planes of corner sharing triangles topologically equivalent to a kagome 
lattice. Powder samples and single crystals with R — Nd and Pr were prepared and 
magnetostatic measurements were performed under magnetic field up to 10 T in the 
temperature range from 1.6 K to 400 K. Analysis of the magnetic susceptibility at 
the high temperatures where only the quadrupolar term of the crystal electric field 
prevails, suggests that the Nd and Pr magnetic moments can be modeled as coplanar 
elliptic rotators perpendicular to the three fold axis of the crystal structure that 
interact antiferromagnetically. Nonetheless, a disordered phase that can be ascribed 
to geometric frustration persists down to the lowest temperature which is about 25 
times smaller than the energy scale for the exchange interactions. 
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The langasite series, the prototype of which is the La3Ga5SiOi4 compound (LGS), 
hence the acronym, belongs to a vast family of materials having the Ca3Ga2Ge40i4 type 
structure [Tj. These materials rapidly attracted a strong interest because they show 
piezoelectric properties with better electro-mechanical coupling and weaker impedance 
than quartz or lithium niobate and tantalate j2] . They can be grown as high quality large 
single crystals, mainly by the Czochralski method, and are now used in surface acoustic 
wave filters in telecommunication devices and high temperature sensors. Crystallizing 
in a non centrosymmetric structure they exhibit quadratic non linear optical and 
electrooptical properties which are also intensively investigated jS] • On the other hand, 
probably because the interest was focused on the striking piezoelectric properties, the 
magnetic behaviors of these materials have not been studied until now, although several 
compounds contain arrays of magnetic cations. We report herein the magnetostatic 
properties of NdaGasSiOu (NGS) and Pr3Ga5SiOi4 (PGS), which are isostructural to 
LGS. Inspection of the crystal structure and analysis of the magnetic susceptibility show 
that these provide us with a rare opportunity for a thorough investigation of geometric 
frustration in kagome magnets, all the more as large single crystals can be grown. 

Geometric frustration is currently attracting strong attention for the numerous 
novel phenomena it might generate |4j. It is classically featured by macroscopic 
degeneracies which prevent magnetic order to set up, thus allowing for new 
states of matter and unconventional excitations to come out 0. Unfortunately 
the materialization of these generally meets with difficulties because of secondary 
interactions (next neighbors or antisymmetric exchange and magnetoelastic interactions) 
or structural and stoichiometric imperfections, which induce ordering at low 
temperatures. Within the temperature range we considered, this is not the case with 
NGS and PGS. 

Langasite type compounds crystallize in the trigonal space group P321, with lattice 
parameters a = b = 8.07 A and c = 5.06 A (for NGS). The structure given in Table^for 
NGS and shown in FigQ contains four crystallographically non-equivalent cation sites, 
and the general formula can be written as A3BC3D2O14. The A site (3e at (~ 0.42 
0)) is 8-coordinated by oxygen anions forming a distorted square antiprism. The B site 
(la at (0 0)) is octahedrally coordinated. The A and B site coordination polyhedra 
share edges to form a layer at z = 0. The C site (3f at (~ 0.76 1/2)) and the D site 
(2d at (1/3 2/3 ~ 0.53)) are both in tetrahedral coordination, the C site being larger 
than the D one. The CO4 and DO4 tetrahedra share corners to form a layer centered 
at z = i, and are also connected to the (A, B) cation layers above and below via corner 
sharing. In the NGS and PGS compounds [2], the rare earth trivalent cations occupy 
the large A sites, while the Ga^+ cations occupy the B, C and half of the D sites. The 
Si^^ cations are localized in the remaining half of the D sites. In Fig|2l the structural 
arrangement of the magnetic rare earth cations in NGS and PGS is outlined. It is 
formed by planes of corner sharing triangles perpendicular to the c-axis and separated 
from each other by the layer of tetrahedral C and D sites. The triangles are equilateral 
by symmetry, being centered by a 3-fold axis passing through the 01 and the mixed D 
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Table 1. Structural parameters for Nd3Ga5SiOi4 determined by single crystal X- 
ray diffraction on a Bruker-Nonius kappaCCD diffractometer (AgKa radiation) on a 
sphere of O.f mm radius, cut from a crystal grown by the floating zone technique. Data 
refined using the SHELX software 6 . Rl= 0.0273, wR2 = 0.0495, GooF = f .059. Cell 
parameters a = 8.066(1) A, c = 5.062 A. The Ga3/Si3 site is half occupied by each 
cation. 
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Gal 


la 











0.01092(9) 
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2/3 
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0.8042(8) 


0.0158(6) 


02 


6g 


0.5341(4) 


0.8514(3) 


0.6916(6) 


0.0201(4) 


03 


6g 


0.2236(4) 


0.0771(4) 


0.7610(5) 


0.0194(4) 



sites, and the intercationic distances are 4.18 A for NGS and 4.21 A for PGS. Exchange 
interactions are mediated by the 01 and 02 oxygen anions. For NGS, the Nd-01-Nd 
angle is 110° with two equal Nd-01 distances at 2.554 A, and the Nd-02-Nd angle is 
106.5° with two Nd-02 distances at 2.38 A and 2.83 A. Because of the lack of hexagonal 
symmetry, this is not an ideal kagome lattice but considering only the shortest atom 
bridging interactions we get the same overall topology. 

Powders of the NGS and PGS compounds have been prepared by solid state 
reactions of stoichiometric amounts of high purity oxides at 1420° C in air. X-ray 
powder diffraction patterns indicated that the samples were single phase and Rietveld 
refinements confirmed the crystal structures reported in Ref. [2]. Single crystals of size 
up to 40 mm in length by 5 mm in diameter were grown by the floating zone method 
using an image furnace, under a 99% Ar + 1% O2 atmosphere, at a growth rate of 10 
mm/h [Zj. Although it does not allow to grow very large crystals, this technique has 
the advantage of preventing a crucible pollution due to long crystal growth times at 
high temperature necessary for the Czochralsky or Brigman techniques. Moreover the 
size of the grown crystal is large enough to allow probing the spatially resolved dynamic 
spin-spin correlations by inelastic triple axis neutron scattering [Hj. The unit cell and 
crystallographic orientations were checked using Laue photographs. The structure of a 
small piece of a NGS crystal was re-determined using single crystal X-ray diffraction 
with a kappaCCD x-ray diffractometer using AgKa radiation. The results are in grood 
agreement with those listed in Table ^ and those reported in j2] . We shall emphasize 
that to the accuracy of the experiment, the Nd site is fully occupied and no substitution 
of another cation for Nd^"^ is detected. 

The magnetostatic properties of NGS and PGS were investigated on single crystals 
under magnetic field up to 10 T in the temperature range from 1.6 K to 400 K on a 
purpose-built magnetometer using the anti-Helmholz two-coil axial extraction method 
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and on a commercial Quantum Design MPMS SQUID magnetometer. We show in Fig|31 
the magnetic isotherms measured on NGS at 1.6 K and 100 K when the magnetic field is 
applied parallel (||) and perpendicular (±) to the three fold axis c of the crystal structure. 
We also give details of the thermal variation of the inverses l/x|| cind l/x± of the 
initial magnetic susceptibilities deduced from the initial slope of the magnetic isotherms 
measured at different temperatures and from the thermal variation of the magnetization 
measured under a magnetic field of 1 T for the same two field orientations with respect 
to the c axis. Similar measurements on PGS are displayed in Fig|3 A large uniaxial 
magnetocrystalline anisotropy is evidenced in both compounds in the whole temperature 
range of the measurements with the c axis as the magnetization axis at low temperature 
changing to a hard axis on increasing the temperature at about 33 K in NGS and 127 
K in PGS. A smaller magnetocrystalline anisotropy in the plane perpendicular to the c 
axis exists which however dwindles out rapidly on increasing the temperature. 

Analysis of the magnetic susceptibility can easily be performed at high temperature 
(T) by limiting the expansion of the quantum statistical average of the component of 
the magnetization along a quantization axis a (M^) to low orders in 1/T as done in jH] 

„ lACi/', 1 rr[i/'0L(/)] , ^,1 

= iv § - lir j( J + 1)(2 J + 1) ' + Q' W 

is the number of the ions in the crystal. C = gjfi'^J{J + l)/3kB is the Curie constant. 
J = 9/2 for Nd^"*" ions and J = 4 for Pr^+ ions, if^ is the component along a of the 
magnetic field on the i-th ion, including molecular field contribution. accounts for 
the crystal electric field potential on the i-th ion and can be expanded as 

= E E i^LYoUJ) (2) 

k q=—k 

in terms of Stevens equivalent operators Ol^{J) [101 115 • A; < 2/ = 6 for f electrons 
(/ = 3). Using the following identity \, 



Tr 



OW)OW) = Tr {3Jl - {Jf f 



= ^ J(J + 1)(2J + 1)(2J - 1)(2J + 3)6k,2Sg,o (3) 

equation is considerably simplified and allows to deduce the inverse of the initial 
uniform magnetic susceptibility Xa along a as 

1^1 (2J-l)(2J + 3) 1 " „ 

where the paramagnetic Neel temperature 6 < accounts for the antiferromagnetic 
exchange interactions. The (^2a)* coefficients are deduced from the crystal symmetry 

J This identity is derived on observing that the product Dk^Vi of two irreducible representations T^k 
and T>i of the rotation group SO (3) contains the trivial representation Vq solely if fc = Z and that the 
trace of Ol^{J)Ol^{J) cancels unless q = — m since under a rotation uj about the a axis this product 
is multiplied by e-'(9+™)'^. 
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at the i-th ion position, which consists of only a two fold axis : = a for the ion at 
(x 0), = h for the ion at (0 x 0) and = u = a + h for the ion at (-x -x 0), 
where x = 0.41809(2) for NGS (see Table P). On choosing as the quantization axis, 
adopting the three fold c axis as the common Y axis for all the ions and completing 
with the appropriate X* axis so that (X\ Y, Z*) forms a right handed frame, V"^ writes 
-^2^2z('^) + -^2^2z('^)+ terms of order 4 and 6 in fc, where B2 and i?| do not depend on 
the i-th ion since, with respect to the (X*,y,Z*) frame, the same crystal environment 
is seen by the i-th. ion. We deduce on rotating the Stevens operators [T2j that 

]^TXAl\{)' = -\iBl + Bl) (5) 
^j:{Ali.r = -ABl + Bl) (6) 

i=l 

We observe first that x± does not depend upon the orientation of a within the 
plane perpendicular to c, which means that the in-plane anisotropy observed at low 
temperature should arise from higher order terms in the crystal electric field potential. 
In NGS, a linear fit of above 100 K and of 1/X|| above 300 K using the slope fitted 
from l/x± yields a paramagnetic Neel temperature 9 = -52 K, an effective moment 
/ieff = 3.77 fiB, close to the value of the Nd^"*" free ion, and a quadrupolar electric 
field parameter {B^ + B2)/kB = -6.35 K. No such quantitative analysis is possible in 
PGS because the temperatures at which l/x\\ and l/x± are linear and parallel to each 
other are beyond the experimental range, but negative 9 and {B2 + -BD/fc^ should be 
expected. As a matter of fact (i?2 + -Bf) < in both NGS and PGS since x\\ < X± at 
high temperature. 

We computed both B2 and -B| for NGS in a point charge model considering the 
eight oxygen anions coordinating the Nd'^"^ ion and forming a distorted square antiprism 
around it. Generally the as-computed absolute values differ significantly from the actual 
ones, but the Bl/B^ ratios are better determined and a much greater confidence should 
be put on them. We get B2/B2 ~ 3/2, which suggests that the Nd magnetic moments 
behave at high temperature most probably as coplanar rotators perpendicular to c with 
an orientational preference along the axis for the i-th ion within the (X*, Y) plane. 
A more complex behavior should occur at lower temperature where higher order terms 
of the crystal electric field compete with the quadrupolar one. 

A fact of outmost importance was that, although finite values of the paramagnetic 
Neel temperature 9 are deduced, no anomaly is detected in the magnetic susceptibilities 
of both NGS and PGS indicating that a disordered phase would persist down to the 
lowest temperature despite antiferromagnetic interactions of one order of magnitude 
larger in energy scale. An elastic neutron scattering test experiment was performed 
using a powder sample of NGS on the D20 high flux diffractometer at the Institut Laiie 
Langevin § that confirms the absence of any magnetic ordering. As evidenced in FigEl 

§ |http:// www, ill, fr / YellowBook /D20| 
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no increase in the Bragg intensity nor emergence of new Bragg intensity are observed 
on the neutron pattern on decreasing the temperature from 100 K down to 2 K. 

In short, geometric frustration in Nd3Ga5SiOi4 (NGS) and PrsGasSiOu (PGS) 
appears to inhibit condensation into a Neel phase and to favor a disordered phase 
down to the lowest temperature. We recently performed inelastic neutron scattering 
experiment on both polycrystals and single crystals evidencing liquid like dynamical 
spin-spin correlations with unusual dispersive features jS] . As from the present study of 
NGS and PGS, we expect that the series of R3Ga5SiOi4 compounds with the other rare 
earth R should show valuable magnetic behaviors inherent to geometric frustration. As 
such they are of as much interest as from their piezoelectric properties. 
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Figure 2. Atomic arrangement of Nd3Ga5SiOi4 projected along the c-axis. The full 
lines linking the Nd^+ cations enhance the magnetic net topologically equivalent to a 
kagome lattice. 




100 200 300 400 

T(K) 

Figure 3. Magnctostatic data measured on a NGS single crystal under magnetic field 
parallel (||) and perpendicular (_L) to the three fold axis c. Top : Magnetic isotherms 
at 1.6 K and 100 K. Bottom : Thermal variation of the inverses l/x|| and 1/X-L 
of the magnetic susceptibilities. The dashed lines stand for high temperature linear 
extrapolations for each field direction. Inset : zoom on the cross-over of l/xy and 
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Figure 4. Magnctostatic data measured on a PGS single crystal under magnetic field 
parallel (||) and perpendicular (_L) to the three fold axis c. Top : Magnetic isotherms 
at 1.6 K and 100 K. Bottom : Thermal variation of the inverses l/x|| and l/x± of the 
magnetic susceptibilities. 
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Figure 5. Neutron patterns collected at 100 K and 2 K from a powder sample of 
NGS and difference pattern multiplied by 10. The consecutive positive and negative 
intensities are due to thermal expansion. 



